function [lon_moor,lat_moor] = Z_tidegaugeLocations(mloc)
%
% A function to store lon, lat for specified tide gauges, and provide that
% information upon user query. Very similar to Z_mooringLocations.m. At 
% present, only tide gauges between latitudes of 43 N and 50 N appear below
% (within the current Cascadia model domain).
%
% The assumption is that tide gauge station positions do not change.  If
% they do and aren't renamed then date information will need to be added 
% (see e.g., Z_mooringLocations.m).
%
% For US NOAA tide gauge data and metadata see:
%     http://tidesandcurrents.noaa.gov/
%
% For Canada tide gauge data and metadata see:
%     http://www.meds-sdmm.dfo-mpo.gc.ca/isdm-gdsi/twl-mne/index-eng.htm
%
%
% INPUT:
%   mloc is a string, e.g. 'La Push' for La Push, WA; 'South Beach' for
%       South Beach, OR; 'Seattle' for Seattle, WA; etc. See below.
%
% OUTPUT:
%   lon_moor [degrees] tide gauge longitude, -180 to 180 format
%   lat_moor [degrees] tide gauge latitude
%
%
% ALTERNATIVELY:
%   A quicker alternative is to simply load the .mat file
%   'tidegaugeLocations.mat' created from 'make_tidegaugeLocations.m' which
%   simply calls this function to make the .mat file.
%
%
% R. McCabe 04/2014


switch mloc
    
    % Salish Sea stations
    case 'Campbell River'
        lat_moor = 50.042; lon_moor = -125.247;
    case 'Point Atkinson'
        lat_moor = 49.337; lon_moor = -123.253;
    case 'Vancouver'
        lat_moor = 49.287; lon_moor = -123.110;
        
    % New Westminster is reasonably far inland up the Fraser River, so ignore...
    %case 'New Westminster'
    %    lat_moor = 49.2; lon_moor = -122.91;
    
    case 'Patricia Bay'
        lat_moor = 48.6536; lon_moor = -123.4515;
    case 'Victoria'
        lat_moor = 48.424; lon_moor = -123.371;
    case 'Cherry Point'
        lat_moor = 48 + 51.8/60; lon_moor = -122 - 45.4/60;
    case 'Friday Harbor'
        lat_moor = 48 + 32.8/60; lon_moor = -123 - 0.6/60;
    case 'Tacoma'
        lat_moor = 47 + 16.0/60; lon_moor = -122 - 24.8/60;
    case 'Seattle'
        lat_moor = 47 + 36.1/60; lon_moor = -122 - 20.3/60;
    case 'Port Townsend'
        lat_moor = 48 + 6.7/60; lon_moor = -122 - 45.4/60;
    case 'Port Angeles'
        lat_moor = 48 + 7.5/60; lon_moor = -123 - 26.4/60;
    case 'Neah Bay'
        lat_moor = 48 + 22.0/60; lon_moor = -124 - 36.7/60;
    
    % Coastal stations
    case 'Tofino'
        lat_moor = 49.154; lon_moor = -125.913;
    case 'Bamfield'
        lat_moor = 48.836; lon_moor = -125.136;
    case 'La Push'
        lat_moor = 47 + 54.8/60; lon_moor = -124 - 38.2/60;
    case 'Westport'
        lat_moor = 46 + 54.2/60; lon_moor = -124 - 6.3/60;
    case 'Toke Point'
        lat_moor = 46 + 42.4/60; lon_moor = -123 - 58.0/60;
    case 'Astoria'
        lat_moor = 46 + 12.4/60; lon_moor = -123 - 46.1/60;
    case 'Garibaldi'
        lat_moor = 45 + 33.2/60; lon_moor = -123 - 55.1/60;
    case 'South Beach'
        lat_moor = 44 + 37.5/60; lon_moor = -124 - 2.5/60;
    case 'Charleston'
        lat_moor = 43 + 20.7/60; lon_moor = -124 - 19.3/60;
        
    %case 'Port Orford'
    %    lat_moor = 42 + 44.3/60; lon_moor = -124 - 29.9/60;
    
    otherwise
        error([mloc ' is not defined in Z_tidegaugeLocations.m'])
end
